Subthreshold behavior and avalanches in an exactly solvable Charge Density Wave system 
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We present a simplified charge density wave (CDW) model exhibiting a depinning transition with 
threshold force and configurations that are explicit. Due to the periodic boundary conditions imposed, 
the threshold configuration has a set of topological defects whose location and number depend on the 
realization of the random phases. These defects are responsible for the randomness of the threshold force, 
which, as the system size L — > cxa, approaches a Gaussian with variance scaling as L^^. The approach to 
threshold proceeds via avalanches whose size dependence on the external driving force F is described by 
a record-breaking process. We find that the depinning transition in this model is a critical phenomenon, 
with the cumulative avalanche size (polarization) diverging near threshold as (Fth — F)~'^ . 
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The motion of Charge Density Waves (CDW) belongs to 
a class of systems in which a deformable structure (elastic 
manifold) is driven by external forces through a medium of 
impurities (random medium). It has been argued that the 
depinning transition for CDWs is a dynamic critical phe- 
nomenon, a phase transition between the pinned and sliding 
states with the driving force as the control parameter [IJ. 
Analytical results for the divergence of strains at the depin- 
ning transition f^, functional renormalization group calcu- 
lations |3], and extensive numerical simulations of CDWs 
and similar systems in 1-3 dimensions fllsl] support this 
claim. However, apart from the functional RG calculations, 
which are perturbative, and simulations, there is little evi- 
dence whether the depinning transition is indeed a critical 
phenomenon, particularly in d = 1. To address this ques- 
tion we introduce an exactly solvable Id CDW model that 
exhibits the tell-tales of a critical phenomenon and allows 
us to understand the origin of criticality. The model is eas- 
ily extended to higher dimensions but we will not pursue 
this in the present letter. 

We begin with the CDW Hamiltonian: 



n{y^}) 



Y.2^y^ 



y.~if + V{y,-a,)-Fy,, (1) 



where V{x) is 1-periodic and a^, the impurity phases, are 
i.i.d. uniform on (— i, +1). Following Narayan and Fisher 
yfl, we choose 

V{x) = {\/2){x-lx\f. (2) 

Here A is the strength of the potential and \x\ is the nearest 
integer to x. Write 

m.i = Ivi - ai\ and y^ = y., - a^ - m„ (3) 
for the well number and well coordinate, respectively, of 
yi, the former records which parabolic well contains yi and 
the latter the displacement G [— 1,+^) from the well's 
center. For large system size L with periodic boundary 
conditions, the static configurations are given as 

^ Y, r?l^-^l (Aa, + Am,) + F/A, (4) 
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where < ?7 < 1 is defined via 2 + A = r/^^ + r], and 
we will ignore the 0{r]^) error introduced by summing in- 
stead over a single period. The nonlinearity of the sys- 
tem is an admissibility condition for the rrii : dU must give 
yi G [— i, +1). All static configurations are linearly stable 
unless a particle is at a cusp of V. 

From (lU we see that if the force F is increased or de- 
creased, the particles will translate uniformly. F can be 
increased until the first particle, say j, reaches a cusp, 
jjj = 1/2. Any further increase in F causes particle j 



to jump wells rrij 



m. 



+ 1, displacing 



yk->yk- Sjk + 



1 






(5) 
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which may cause other particles to jump also. We assume 
that the dynamics of the particles are purely relaxational 
and that all changes in the force are sufficiently slow that 
we reach the static configurations. 

Like the CDW models with sinusoidal potential [51], this 
model exhibits reversible and irreversible behavior. If we 
change the force by AF, causing one or more particles to 
jump wells before reaching another stable configuration, 
and then return F to its original value, it can be shown la] 
that none of the particles jump back. On the other hand, 
the resulting configuration moves uniformly and without 
jumps, i.e. reversibly, in the interval of forces between F 
and F + /S.F. This establishes an equivalence relation 
among configurations that share the same well numbers ttt-j 
but differ in their well-coordinates y^ by an overall transla- 
tion. We choose as representatives the zero-force configu- 
rations satisfying (|4]l with F = 0. 

We simplify further by choosing r] small: defining the 
rescaled well coordinate Zi as rjZi = yi — F/X, we obtain 
fromdU 



Zt 



Aa, + Am, + 0(1]) 



(6) 
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We will refer to the CDW model, with 0{rf) terms 
dropped {0{r]) in the rescaled coordinates), as the toy 
model. A set of exact results can be obtained for the toy 



model; the proofs and additional results for the untruncated 
model will be given elsewhere Ig]. 

It is readily shown that the process of increasing the ex- 
ternal force and evolving the configurations to threshold 
can be described by the zero-force algorithm (ZFA), which 
iterates the following procedure: 

(i) Record z^ax = max, z,. 

(ii) Find j = argmax^Zi JtIi, set Zj -^ Zj — 2, 

Zj±i —7- Zj±i + 1, and rUj — )• rrij + 1. 

(iii) If any z, > z^ax from (i), goto (ii). 

The first execution of (ii) jumps the particle which would 
first reach the cusp if the force were increased (uniformly 
translating the particles) and subsequent executions resolve 
those particles which would be pulled across the cusp as a 
result of the first. This process always terminates, possibly 
with all sites having jumped exactly once, in which case 
the Zi are unchanged. This fixed-point is the threshold con- 
figuration, and the ZFA finds it in finitely many steps. 

From the ZFA it is clear that the evolution to the depin- 
ning transition minimizes niax^ Zi and indeed the thresh- 
old configuration is the solution of this optimization prob- 
lem; this is true for the non-truncated CDW model as 
well. Equation ^ suggests that setting Ami = ~l^(^i} 
would be favorable, but periodicity and the requirement 
that rrii G Z usually prevent this. The threshold config- 



L-lr 



uration depends on S* = J2i=o l^'^t} ^^'^ i^ given as 



Am, = -lAa,j + J, 



JikL 



(7) 



where J is an integer vector selected as follows: 



(i) Case S" > 0. Ji = 1 for the S" + 1 positions i which 
have smallest Aa^ — JAaJ, J^ = otherwise; 

(ii) Case S < 0. J^ = — 1 for the l^l — 1 positions i 
which have largest Aa^ — [AaJ, J, = otherwise; 



and ki, is an index defined by (divisibility condition) 

-lAa,i + Ji) (modL). (8) 



L-l 

4=0 



We will refer to the sites where Am,i + \Aai 
as defects with charge e^. 

Given the threshold configuration zf^, the corresponding 
threshold force is 



i^th ({a.}) = A (1/2 - ?7z: 



th \ 
max/ ' 



(9) 



where z^^x = max^ zf^. The term in parentheses on the 
RHS is the distance from the particle with maximum z^ to 
the cusp. Writing U; = Aa^ — [AaJ and using (|6]l and 
(|7]l, we find 






(10) 
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FIG. 1. Evolution of the m^ from negative (blue) to positive 
threshold configuration (top red curve) for L — 512. Approach- 
ing the threshold, increasingly large segments depin and advance. 
The inset shows the last configuration change Srrii — m^ — mi 
and the corresponding labeling of the avalanche series. Note the 
trapezoidal shape of Srrii: the comers satisfy ASnii — ±1. 



where a is the permutation of the indices that orders uj: 



w, 



^(0) 



< ^r 



(1) 



< 



< W^(L-1 



The results for the 



threshold configuration and force characterize the manner 
of their dependence on the random phases q^, via (a) the 
statistics of S and (b) the rank statistics of {cOi} [8]. Rou- 
tine arguments show L^^^^S converges in distribution as 
L — 7- oo to a normal random variable with mean and 
variance 1/12, so the typical number of topological de- 
fects scales as L^^^. As a consequence the variance of 
the threshold force distribution averaged over the quenched 
disorder is found to scale as L^^ and thus fluctuations 
AFth around the mean EFth = (1 - r]f/2r] + 0{L-^''^) 
scale as L~^/^. This matches the scaling behavior of the si- 
nusoidal CDW model |15|], for which the expected behavior 
is AFth ~ L^^/"^, with Vt the finite-size scaling expo- 
nent from the scaling theory of Chayes ef a/ [|9|]. Our model 
saturates their prediction i^^ > 2 in Id. 

Equation (|7]) can be used to show that if the strains Si = 
rriij^i — rrii are rescaled as 



\l/2 



s^^\t) ^ {12/ Lf'\sy,/^^ (0<t<l), 



(11) 



the resulting cadlag processes s^^^ converge in distribution 
as L — > oo to a periodic Brownian motion with zero inte- 
gral: 



B{t) 



B{r)dr (0 < t < 1), 



(12) 



where B{t) is a standard Brownian bridge. It follows that 
the typical maximum and minimum of (sj) scale like L^/^, 
matching our expectation for diverging strains at threshold 
as shown by Coppersmith [2]. 

We now consider the behavior of the system as it ap- 
proaches threshold. As in ^ we start out with the con- 
figuration {Zj~} at negative threshold force, F^^ < 0, and 
increase the force until we reach the positive threshold con- 
figuration {Zj^}. Figure [T] shows the corresponding evo- 
lution in terms of the well numbers {rrii}. The negative 
threshold configuration {z,^} can be obtained from (|7]l by 



switching a — )• — a. For simplicity we consider the case 
5 = 1 for which 

z'l = uj., + (5j,fe-, (13) 

Z+ =UJi + 6i^a(0) + ^i,a{l) - SiM , (14) 

with k"^ determined by ^. Following the ZFA, the first 
site to jump will be k^ and z,nax = ^k- + 1. As a result 



Zk- 



bJk- + 1 



^k- 



1 



and 



Zk-±l = Wfc-±l -^ Wfc-ii + 1. 



The neighboring sites A;~ ± 1 will be forced to jump if 
ujk-±i > ujk-- (Here and in the following addition and 
subtraction of indices are mod L.) It is readily shown that 
all those consecutive neighbors j to the left and right of 
k~ for which uJj > uj^- are forced to jump, giving rise 
to an avalanche with extents i^ and ij^, respectively, where 
Jl = k" —iL — lis the first site to the left of k^ for which 
LOj^ < UJk- and likewise for jj^ = k^ + II + 1- Note that 
if a site j and both its neighbors jump, Zj is unchanged; in 
particular if A:~ it 1 jump, then Zk- is still the maximum 
and will initiate the next avalanche. 

The above observations let us group several iterations in 
the ZFA which resolve all the avalanches initiated at a sin- 
gle site. To understand this process, observe that only the 
ranks a{j) (not the precise values uJj) are required to deter- 
mine the avalanche waves. From Zi we draw a rank-defect 
diagram, displaying the rank a{i) of oj{i) with over- or 
underlines to indicate positive or negative defects, respec- 
tively. Consider an example: let z^ have representation 

... 10 12 17 15 16 18 13 14 1 .. . 

so that cr(k^) = 15. The extent of the first avalanche 
will he i^ = l,ifl = 2 and after the bracketed sites have 
jumped, the resulting configuration becomes 

. . . 10 T2 [17 T5 16 18] TT 13 1 . . . . 

Next k^ and k^ + 1 will jump again, yielding 

... 10 T2 17 [15 16] 18 TT 13 1 . . . , 

No more avalanche waves are initiated at rank 15. Gener- 
alizing, given an initial site j, and initial extents {iL^iii), 
a new configuration will be reached after min(ii, i^) + 1 
avalanches, at the end of which the sites at j^ , Jb. will have 
acquired a charge +1, while the sites k^ and k^ + i^ — i^ 
will each acquire a charge —1. As a result, a total number 
of (ii + 1) {in + 1) jumps will have occurred. The next se- 
ries of avalanches will be initiated by the rank 12 site with 
initial extent (0, 2) and the configuration 

... To [12 17 15] 16 18 TT 13 1 .. . 

is reached after 1x3 = 3 jumps. Next, the rank 11 and 10 
sites trigger avalanche waves of extents (2, 1) and (0, 4), 
respectively, and we finally reach 

... 10 12 17 15 16 18 11 13 T .... 
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FIG. 2. Numerical results for the expectation of the cumulative 
jump size (polarization) plotted against the reduced force. Plot 
symbol colors refer to different system sizes as indicated in the 
legend. The numbers of realizations are indicated in parenthe- 
ses. Scaling of the axes as indicated in the figure. The horizontal 
dashed line corresponds to dlSt , the slanted dashed line indicates 
a power law of exponent ~2. The solid line is the finite size scal- 
ing function i20\i . 



This is the threshold configuration as predicted by (IT4l : 
the rank 1 particle would jump next, and the extent of this 
avalanche would be (3, L — 5) (periodic BC). The rank 16 
site would have both neighbors jump and acquire a net pos- 
itive defect, giving it a value above z^ax = '^^(i) + 1; and 
it must also jump. Thus all particles jump, and the config- 
uration is unchanged. 

We now consider the statistics associated with this 
threshold-to-threshold evolution. Given a sequence of ran- 
dom variables Xi, X2, . . ., the site i is said to be a lower 
record, if min(Xi,X2, . . . ,Xi) = Xi. As indicated in the 
example, the sites initiating the avalanches are the lower 
records of the sequences 

and Jii = u)k-,^k-+i,^k-+2-,---,^aR, 

where {o"/^, CTt^} = {o"(0), o"(l)} are the termination sites 
iflOll . The evolution from negative to positive threshold ter- 
minates when avalanches initiated dA k^ and the left and 
right lower records reach a^ and Gr. 

Observe that the active region where particle jumps oc- 
cur is given by the segment containing k^ with boundaries 
Gl and Gj^, and that the number of jumps leading to the rth 
configuration, whose left and right boundary have reached 
3l{t) > Gl and Jr{t) < Gr, is 



S(r) 



E 



m,{T) 



m. 



{k- - uir)) (Mr) - k-) . 



(16) 



In particular the total number of jumps in the threshold-to- 
threshold evolution is 

j: = {k- - gl) {(Tr - k-) . (17) 

The random variables oji are exchangeable and statistics 
of records are well known in this case [11]. From this the 



statistics of S and E^- can be obtained. We defer details to 
la] and focus here on expected values. By exchangeability 
a"i and 0"^ are selected uniformly from all (distinct) pairs 
of indices, and [6] k~ is uniform and independent of a. 

The expected total number of configurations N^^^p^ in 
the threshold-to-threshold evolution is related to the ex- 
pected number of records in the sequences (ITSt , and one 
finds EA'stcps ~ 2 In L — 4 + 27, where 7 is Euler's con- 
stant. The distribution of the total number of jumps S can 
be obtained from (ITtI i and the statistics of k^, Gl and Gr 
^. In particular, 



ES = LV12 + 0(L). 



(18) 



Numerical results for the mean of the polarization (IT6] 



as a function of the reduced force F^ 



,th 



Fl ^ z^ 



.(r) 



, cf. dUl, are shown in the scaling plot Fig. |2l The 
data collapse for different sizes L is perfect. To explain the 
scaling behavior and cross-over between the two regimes, 
we first require some notation. 



Write X 



^th 



€ [0, 1) for the maximum 
height of the current configuration less that of the (+)- 
threshold configuration. Shift indices so that k~ = 0, 
and let Jl{x) and Jb{x) be the (noninclusive) left and 
right extents of the interval of sites that have jumped in the 
threshold-to-threshold evolution to achieve X < x, cho- 
sen so that —L<Jl < < jfl < Jl + L. Then the 
corresponding polarization is S^ = —Jl{x)Jb{x)^ where 
Jl (x) and j^. (x) are the first indices in JTl and J'ji of (ITSll. 
respectively, for which uJi < <^ct(i) + x. 

For large L, approximate (^^(i) = —\ and w in J^ and 
Jr as Ltd. uniform (— 1,+|) variates |0], sharing their 
first elements. Their lengths 1.^ and ^^ are random, inde- 
pendent of the sequence values, and uniform over all possi- 
ble values. The random variables \Jl \ and j^ are geomet- 
ric, truncated at II and (.r, and an easy calculation gives 
^-3l{x)Jb{x) KL,£fl]: 

1 



"(1 



1 



•)(1 



1 



(19) 



Setting X = u/L and averaging over i^ and i^, 
we obtain the finite-size scaling function $(t() = 

limi^^L-2E[S„/i]: 



$(n) 



6 - 4n + n^ 



6e" 



2ne"" 



n-* 



(20) 



This has no pole at u = 0; expanding the exponentials, 

$(u) = l/12-u/30 + u7l20-uVl05 + O(u'^) (21) 

for < u <C 1, and for n ^ 1 we have ^[u) ~ n^^, 
with a crossover near u = 1. More intuitively, Jb{u/L) is 
a geometric random variable with mean L/u — 1 truncated 
at £fl, which has mean L/3 + 0(1). When u ^ 1 the 
truncation is largely irrelevant, since Jr <C Ir with high 
probability, and Jr{u/L) scales as u~^. As u nears 1 the 
truncation can no longer be ignored, and it dominates for 
u <^ 1, where Jr{u/L) is insensitive to the exact value 



of u. In terms of finite-size scaling, for each of —j^ {x) 
and Jr{x) we have a correlation length ^ ~ x^" with 
finite-size scaling exponent v = 1 and finite-size effects 
dominating when L/^ <^ 1. Even though j^ and Jr are 
dependent, the dependence is not so strong as to impact the 
scaling exponents of the product: they still combine to give 
scaling exponent —7 + 1 = —2 as in Fig.|2] 

We have presented an exactly solvable Id CDW toy 
model with a critical depinning transition, exposing the 
role played by the disorder and the boundary conditions. 
The approach to threshold is governed by a stochastic 
record-breaking process. The scaling exponents of the toy 
model differ from those obtained in terms of the 4 — e ex- 
pansion of Narayan etal. [IsId], in which u = 2 and 7 = 4 
are found in d = 1, but agree with those of the d = 1 CDW 
automaton model of Myers and Sethna [4]. These authors 
also note that for d > 1 the automaton model yields ex- 
ponents in agreement with those obtained from the 4 — e 
expansion. We have also carried out numerical simulations 
for the full CDW model, dU, and find that the exponents i^ 
and 7 agree with those predicted by Narayan et al. Hence 
the full and toy CDW model exhibit different subthreshold 
critical behavior. We will report detailed results on the full 
model elsewhere. Lastly, it should be interesting to study 
the scaling behavior of the toy model presented here for 
d>l. 

DK and MM would like to thank F. Rezakhanlou for 
stimulating discussions, MM acknowledges discussions 
with H.J. Jensen, partial support by Bogazi9i University 
grant 12B03P4 and thanks the Berkeley Math department 
for their kind hospitality during his sabbatical stay. 



t 



kaspar® math.berkeley.edu | 



mmungan@boun.edu.tr 
[1] D.S. Fisher, Phys. Rev. Lett 50, 1486 (1983)Phys. Rev. B 

31, 1396 (1985) 
[2] S. N. Coppersmith, Phys. Rev. Lett. 65, 1044 (1990)S. N. 

Coppersmith and A. J. Millis, Phys. Rev. B 44, 7799 (1991) 
[3] O. Narayan and D.S. Fisher, Phys. Rev. B 46, 11520 

(1992)R Le Doussal, K.J. Wiese, and P Chauve, 66, 174201 

(2002)D. Ertas and M. Kardar, 53, 3520 (1996) 
[4] A. Erzan, E. Veermans, R. Heijungs, and L. Pietronero, 

Phys. Rev. B 41, 11522 (1990)C.R. Myers and J.P Sethna, 

47, 11171 (1993)A. Rosso and W. Krauth, Phys. Rev. E 65, 

025101 (2002)H.J. Jensen, J. Phys. A 28, 1861 (1995) 
[5] A.A. Middleton and D.S. Fisher, Phys. Rev. B 47, 3530 

(1993)0. Narayan and A.A. Middleton, 49, 244 (1994) 
[6] D.C. Kaspar and M. Mungan, in preparation 
[7] We assume the number particles L is large but finite, so that 

such a particle can always be identified. 
[8] Any L — 1 of the {uji} (but not all L) turn out to be i.i.d. 

uniform (-1/2,1/2). 
[9] J.T. Chayes, L. Chayes, D.S. Fisher, and T. Spencer, Phys. 

Rev. Lett. 57, 2999(1986) 
[10] For arbitrary S, a^ and aR are the two sites of the lower 



threshold configuration with lowest z value. 



[11] N. Click, AMM 85, 2-26 (1978)N. Balakrishnan 
B.C. Arnold and H.N. Nagaraja, Records (Wiley, 1998) 



